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Abstract: The chemical potential of the electron gas on a one-dimensional lattice is determined within the 
discrete Hubbard model. The result will have applications in theoretical studies of transport properties of 
quasi one-dimensional organic conductors such as the Bechgaard salts. 
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Introduction 



Quasi one-dimensional (Q1D) organic metals were discovered in 1980 [1], [2], [2a]. It was soon shown 
[3] that the electrical conductivity of these materials could not be described by the theory of electrical 
conductivity of normal metals. The Q1D organic metals are, (since [4]), considered within the framework 
of the Hubbard model on a one-dimensional lattice. In spite of the fact that this model is one of the 
simplest in statistical mechanics, it is a subject of intensive theoretical studies. The model is solvable in one 
dimension, and the chemical potential is known only for the special case of a one-dimensional half-filled 
band, for which \i = [5] .It is usually assumed that the electrons in Q1D organic metals form a normal Fermi 
liquid (however, see [6] ), which means that 

n k = 1 / [l + exp/3(e k - ft)] (1) 

noindent The electronic energies are given by £k = — 2tcos(ks) [5]. All the symbols have their usual 
meaning, t denotes the hopping and s the lattice constant. 

The purpose of the present letter is to determine the dependence on the band-filling and 
temperature of the chemical potential of the electron gas on a one-dimensional lattice. A similar problem has 
been investigated in quantum field theory, where a special prescription had to be proposed for introducing 
the chemical potential in a lattice regulated field theory (see, for example, [15]). The knowledge of the function 
/i = /x(n, /?, t) is of paramount importance for the calculation of the electrical conductivity and optical 
response functions of correlated electron systems [13], [14]. Deviating the filling from 1/2 is a theoretical 
procedure for representing the effects of doping on the material. Preliminary results of this calculation were 
recently published [5a]. 

The calculation which will be reported in this letter is a first step towards a theoretical study of the 
effects of doping on the transport properties and optical response functions of Q1D organic metals. For a 
recent example of an experimental study of these effects see [7]. Apart from the Q1D organic metals, the 
calculation to be reported is of interest for general one-dimensional models of solids ( such as the Kronig- 
Penney model) [8]. 



The calculation 

The starting point for the determination of the chemical potential is the following relation (for example, 

[9]): 

n = 2 J 1/ [1 + cxp (3(e k - /*)] dp / (27rh) (2) 

where all the symbols have their usual meaning, p = hk, and n denotes the number of electrons per site (n 
= 1 corresponds to filling 1/2). The limits of integration are ±7rh/s. Using the fact that in the ID Hubbard 
model £k = — 2tcos(ks), it follows 

dp = h dk = (h/s) (4t 2 - e 2 )- 1/2 de (3) 
Inserting (3) into (2) one gets, after some algebra, that 

2t 

n = 1/(tts) J 1/ [1 + exp /3(e k - ^)](4t 2 - e 2 )- 1 ' 2 de (4) 

-2t 

and the problem is to find a solution of this integral equation. Instead of using the full rigour of the theory 
of integral equations (such as [10]), the problem can be solved by a suitable developement of the Fermi 
function. 

It has recently been shown [11] that the Fermi function can be represented in the following form: 
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1 / [1 + cxp /3(e k - »)] = Q(jm - e) - ^ A 2k+1 /r 2 ( k+1 )<5( 2k+1 ) (e - M ) (5) 



k=0 



where A 2k+1 = (-l) k + 2 [2(2 2k + 1 - l)7r 2k + 2 / (2k + 2)\]B 2k+2 

In eq. (5) 6 denotes the step function, <5( 2k+1 ) (e — fi) are the derivatives of the 8 function and B 2k+2 
are Bernoulli's numbers [12]. Inserting (5) into (4) it follows that 
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n = 1/M [ [9( M -e) - ]T A 2k+1 /J- 2 ( k + 1 )^ 2k + 1 )( £ - M )](4t 2 - e 2 )^ 2 de (6) 



-2t k=0 



This equation can be integrated by using the following relation, valid for an arbitrary function f(x) and 
its n-th order derivative: 



(7) 



yV n )(x-X )f(x)dx = (-l) n f( n )( X0 ) 

Combining (7) and (6), and limiting (5) to terms with k < 2 one finally obtains 
1 / (7rs)[7r + (^ 2 /6/3 2 )a/9e((2t) 2 - s 2 y 1/2 + 7/360(^//3) 4 

d 3 /ds 3 ((2t) 2 -e 2 y 1/2 + 31/(21 x 720)(7r//3) 6 a 5 /9e 5 ((2t) 2 -e 2 )- 1 / 2 ] = n (8) 

where e = \i because of eq.(7). 

Expression (8) can be developed up to the second order in /j, around the point [i = 0, and then solved 
with the following result: 

fj, = (/3t) 6 (ns - l)|t| / [1.1029 + .1694(/3t) 2 + .0654(/3t) 4 ] (9) 
In the case of weak hopping, this can be developed as 

(j, = (/3t) 6 (ns - l)|t| [.90671 - .13940(/3t) 2 - .03284(/3t) 4 ] (9a) 



It is obvious that lim fj, = 0, which means that equations (9) and (9a) reduce to the solution discussed in [5]. 

n,s— >1 

A better approximation of the function [i = t, n) can be obtained by pushing the development of eq. (8) 
around the point fi = to higher orders. As a test, this development was performed to third order. However, 
solutions obtained in this way are too complicated for being applicable in calculations of the conductivity of 
Q1D organic conductors. 



Discussion 



In this letter we have determined the dependence of the chemical potential of the electron gas on a 
one-dimensional lattice on the band-filling, hopping and temperature. The calculation was performed 
within the discrete Hubbard model. It was motivated by the need of introducing the doping in theoretical 
studies of the electrical conductivity of Q1D organic metals. 

The developement of eq. (8) in fj, around the point /x = limits the applicability of eq. (9) to values 
of n not too different from 1. Work is already in progress aiming at an expression for the chemical 
potential which could be applied at arbitrary filling. Another line of progress concerns the application of 
the results of this letter to the analysis of the dependence of the conductivity of Q1D organic conductors 
on the doping. 
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